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Abstract 



The neutron's electric form factor contains vital information on nucleon 
structure, but its interpretation within many models has been obscured by 
relativistic effects. I demonstrate that, to leading order in the relativistic ex- 
pansion of a constituent quark model, the Foldy term cancels exactly against 
a contribution to the Dirac form factor F\ to leave intact the naive inter- 
pretation of G% as arising from the neutron's rest frame charge distribution. 
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I. INTRODUCTION 



In 1962, Sachs showed [1] that the combinations of elastic nucleon form factors (N = p 
or n) 



G« M = F~ + F~ (2) 

have simple interpretations as the spatial Fourier transforms of the nucleons' charge and 
magnetization distributions in the Breit frame (where momentum p — —® is scattered to 
momentum p' = +%)■ Here and are the Dirac and Pauli form factors, respectively, 
defined by 

(N(p', s'M m mN(P, s)) = u(p', S ')[FiV + iTT^FfW, s) (3) 

where q u = pj — p u and Ff and F^ are functions of Q 2 = —q 2 . 

These form factors obviously contain vital information on the internal composition of 
the nucleons. Although it has proven elusive experimentally, the electric form factor of the 
neutron G% is particularly fascinating in this respect. In pion-nucleon theory, G\ would 
arise from a ir~ cloud with convection currents producing the anomalous magnetic moments 
F% = 1.79 = fi p — 1 and i 7 ^ = —1.91 = fi n . In contrast, in a valence quark model, the 
nucleon magnetic moments arise from the underlying charged spin-| constituents with the 
famous 577(6) relation 

= - - (A) 

2 y > 



and with a scale set by 



Hp = — ^ 3 (5) 
m d 



where ~ m u ~ |mjv is a valence quark effective mass. Within this model it has been 
argued that in the SU(Q) limit G^Q 2 ) would be identically zero, but that the spin-spin 



forces which produce the S't/(6)-breaking A — N splitting create a charge segregation inside 
the neutron and lead to a nonzero G E []2|-|J. The effect arises because the spin-spin forces 
push d quarks to the periphery of the neutron and pull the u quark to the center. Thus both 
the 7T~ cloud picture and the hyperfine-perturbed quark model predict a negative neutron 
charge radius, as observed. 

Nonrelativistically, the squared charge radius is simply the charge-weighted mean square 
position of the constituents. More generally 

G P E (Q 2 ) = 1 " \rl P Q 2 + ■■■ (6) 

and 

G n E (Q 2 ) = ~\r\ n Q 2 +■■ (7) 
define the proton and neutron charge radii, with 

= 1 " irV? + " " " (8) 

defining the corresponding magnetic radii. Experimentally, the three form factors G E , G P M , 
and Gm are reasonably well known, although new measurements should soon determine 
them with much greater precision ||. From these measurements, which cover Q 2 in the 
multi-GeV 2 range, we have learned that all three form factors have roughly similar shapes: 

G P E (Q2) ~ G M(Q 2 ) „ G M(Q 2 ) „ (Q) 
ftp ftn 

where the dipole form factor 

°° {Q2) - FTW^P (10) 

with Mj ipole = 0.71 GeV 2 . I note in passing that this observation makes little sense in a 
picture where the nucleons have a point-like core surrounded by a pion cloud (since in such 
a picture G 7 ^ is purely pionic while G P E and G P M are mixtures of a "bare" proton and a pion 
cloud). In contrast, Eq. (9) is very natural in a valence quark model where, to leading order 



in a relativistic expansion, Gd{Q 2 ) would simply be the Fourier transform of the square of 
the ground state spatial wavefunction. 

Recent advances in experimental technique should lead to a clear measurement of the 
neutron's electric form factor G E (Q 2 ) for Q 2 in the GeV 2 range in the next few years 
|J. Its charge radius r 2 En is known from low energy neutron-electron elastic scattering 
to be —0.113 ± 0.005 fm 2 0, but the electron-neutron scattering measurements needed 
to determine G E in the GeV 2 range (in order to roughly map out the neutron's electric 
structure with a resolution of 10% of the proton's size) have been plagued by the lack of a 
free neutron target ||. Fortunately, the method of Arnold, Carlson, and Gross [|J, which 
uses spin observables sensitive to G^ — G^ interference, opens up new methods for measuring 
G E , and recent advances in accelerator, target, and detector technology are beginning to 
exploit these new methods. 

In anticipation of these measurements, there has been renewed discussion about their 
interpretation. I focus here on the belief that the measured r 2 En is explained by the "Foldy 



term" pi. I.e., using Eqs. (1) and (7), 



2 _ 2 3/l n 2 , 2 | | 

' En — ' In ^ ~ ' In ^ ' Foldy, n ) 1 1 1 1 



where r\ n is the "charge radius" associated with F" ~ ~\ r \ n Q 2 + " ' "■ The second term in 
Eq. (11), called the Foldy term, appears to arise as a relativistic correction associated with 
the neutron's magnetic moment and so to have nothing to do with the neutron's rest frame 
charge distribution. It has the value —0.126 fm 2 , nearly coinciding with the measured value. 
On this basis it has been argued that any "true" charge distribution effect must be very 
small. In this paper I will show that while the Foldy term closely resembles r 2 En numerically, 
it does not "explain it". Indeed, I will demonstrate that, in the relativistic approximation 
to the constituent quark model in which the Foldy term first appears, it is cancelled exactly 
by a contribution to the Dirac form factor F\ leaving r 2 En correctly interpreted as arising 
entirely from the rest frame internal charge distribution of the neutron. 
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II. THE INTERPRETATION OF THE NEUTRON CHARGE RADIUS IN A 

CONSTITUENT QUARK MODEL 



The relationship (1) between the Sachs form factor Ge and the Dirac and Pauli form 
factors F\ and F2 is relativistic in origin. Unfortunately, relativistic constituent models of 
the nucleon are notoriously difficult: rest frame models are difficult to boost and infinite- 
momentum-frame (or light-cone) quark models have trouble constructing states of definite 
J p . This could be the reason that the interpretation of G E has not been clarified in the 
context of such models. 

While an accurate constituent quark model of nucleon structure must certainly be fully 
relativistic, the issue at hand can be resolved by using a relativistic expansion around the 
nonrelativistic limit. This is possible because the Foldy term r Fold arises at order Q 2 /m 2 
and so its character may be exposed by an expansion of G E to order 1/m 2 . I will also 
exploit symmetries of the problem available in certain limits which will make the discussion 
independent of the details of models. 

I begin with a simple "toy model" in which a "toy neutron" n§ D is composed of a 
scalar antiquark S of mass and charge —en and a spin-| Dirac particle D of mass 
Trip and charge en bound by flavor and momentum independent forces into a rest frame 
nonrelativistic S- wave. The calculation begins by noting that, from their definitions, 

GT(Q 2 ) = (nsD(+^,+)\Pem\ns D (-^,+)) (12) 

and 

GT(Q 2 ) = ^(ns D (+^ ,+M£ 2 \nsn(-^,-)) • (13) 

It is immediately clear that the calculation of these form factors requires boosting the rest 
frame S-wave bound state to momenta Doing so can introduce a host of 1/m 2 effects 

in the boosted counterpart of the S'-wave state and it can also produce new P-wave-like 
components by Wigner-rotation of the .D-quark spinors |J. I will show that the latter 
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effect is subleading, and will deal with the former effect by exploiting an effective charge- 
conjugation symmetry of the system for m D = m s = m. 

Since Ph §D involves the limit of Eq. (13) as Q — > 0, to the required order in 1/m it 
simply takes on its nonrelativistic value 

where of course mn SD = ms + mp in this limit. The Foldy term is thus well-defined: 

r Foldy,hg D ~ • (15) 

We next compute G n ^ D (Q 2 ) directly from Eq. (12). To leading order in 1/m 2 , p em remains 
a one-body current, and the impulse approximation is valid. Within this approximation, we 
make use of the relation 

Q 2 

(D(p + Qz, s')\p em \D(p, s)) = e D (l - —^ r )(D(p + Qz)\p em \D(p))5 ss , + p sp in-fu P (16) 
where 

Pspin-fiip = -^r(P- s s>+5s- - p+5 S '~5 s+ ) (17) 

and where D is a fictitious scalar quark with the mass and charge of D. This expression is 
easily obtained by making a nonrelativistic expansion of both the D and D charge density 
matrix elements. 

The spin-flip term p sp in-fiip can only contribute via transitions to and from the Wigner- 
rotated components of the wavefunction. However, the amplitudes of such components are 
proportional to Qk/m 2 D , where k is an internal momentum. Since p sp in-fiip already carries a 
factor 1/m 2 }, such effects may be discarded. Note that non-flip Wigner-rotated contributions 
are of the same order and may also be neglected. 

We conclude that r 2 Efl _ D may be computed by replacing D by D provided the additional 
contribution —eBQ 2 j9>m 2 D is added to G^f D . (More precisely, this factor multiplies the 
D contribution to G E SD , but at Q 2 = this is just unity.) I will denote the associated 
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"zwitterbewegung" charge radius 3e£>/4m|, by r 2 D zwitter . The effect of r 2 D zwitter is well- 
known in a variety of contexts, including atomic physics JTTJ, nuclear physics hadronic 
physics ||13|| , and heavy quark physics ||14| ; at the most elementary and concrete level it 



appears as the additional factor of (1 — i^r) 2 in the ratio of the Mott cross section to 
the Rutherford cross section. The problem of computing the remaining contributions to 
G^f D to this order from the fictitious SD scalar-scalar bound state would in general be 
highly nontrivial. However, in the limit tud = m>s = m these contributions vanish, since 
this system has in this limit a pseudo-charge-conjugation invariance under (S, D) — > (S, D). 
However, we note that in this limit m^ D = 2m and so from Eq. (15) we have 

r 2 =r 2 =^JL =r 2 (]Q\ 

' En SD ' D, zwitter ^ ffl 2 Foldy,hg D ) \ ±U J 

i.e., in this model the "scalar charge distribution" is zero and the Foldy term would indeed 
account for the full charge radius of h. This conclusion is simply interpreted: the two scalar 
particles 5* and D have perfectly overlapping and cancelling charge distributions, but the 
expansion of the D distribution by r 2 D zwitter creates a slight excess of D at large radii. 
In terms of its experimental significance, we have concluded that in an SD model of the 
neutron, the observation of an equality of r 2 En and r 2 Foldy n would indeed indicate the absence 
of an intrinsic "scalar" charge distribution. 

We shall soon be drawing quite another conclusion for the situation in the constituent 
quark model. However, before leaving the SD model, it is useful to consider another limit: 
the "hydrogenic limit" where m s — > oo. In this case (see Eq. (15)), r 2 Foldy fl _ D = but 

2 ^ e D ,2 /in\ 

rE ' h sD = + e ^ wf , (19) 

where r 2 w j is the charge radius associated with the bound state D problem. It may be 
that r^j contains other l/m 2 D effects, but we note that the physics of the r 2 D zwitter effect 
previously associated with F 2 via r 2 Foldy fl _ D now must be asociated with F% via r\n §D - The 
SD toy model thus simultaneously supplies us with a simple interpretation of the Foldy 
term and a warning about associating F\ with the neutron's "intrinsic charge distribution" . 



While the SD toy model has some of the characteristics of diquark models for the nucleoli 
(S has the quantum numbers and color of a scalar diquark), our main use for it was to 
introduce the basic elements of our discussion in a simple context. Indeed, it is now relatively 
trivial to extend our considerations to the realistic case of the valence quark model in which 
the neutron is in the leading approximation made of three mass m q spin-| quarks ddu bound 
by flavor and momentum independent forces into flavor-independent nonrelativistic relative 
S'-waves. In this case 

f, = _^^^_2 (20) 
Pdi " 3m q V ; 



so that 



r Foldy,n dd „. ~ • ( 21 ) 



In calculating r% h via Eq. (12), the transformation of the calculation of the charge radius 
of ddu to that of three scalar quarks ddu and residual r 2 zwitter = —e q Q 2 /8m 2 terms proceeds 
as before, as does the neglect of Wigner-rotated components of the boosted state vectors. 
However, in this case, since the r 2 zwitter terms are spin and flavor independent, and since the 
sum of the three charges is zero, they lead to no net Q 2 /m 2 term! The reason for this is clear: 
the exactly overlapping and cancelling quark distributions remain exactly overlapping and 



cancelling after they are all equally smeared by r 2 zwitter Ul5| . This picture also anticipates 



the next stage of the argument: the analog of the pseudo-charge-conjugation invariance that 
we used for the scalar part of the SD matrix element is that the three quark wavefunction 
belongs to the symmetric representation of the permutation group S3 so that the scalar part 
r En--- °f the charge radius vanishes. Thus in the usual valence quark model 



r E,n ddu — r I,n d - d - fi + ^i e i r i,zwitter ~ > ( 22 ) 



which requires that 



r tn ddu r Foldy,n ddu 

(23) 



S 



so it is appropriate to interpret the observed r 2 En as due to an intrinsic internal charge 
distribution. Stated in another way, the coincidence of the predicted rest frame charge 
distribution in such models with the experimental value of r\ n may be claimed as a success, 
while the numerical coincidence of r 2 En with the Foldy term may consistently be viewed 
as a potentially misleading accident. Such an accident is possible because while in the 
nonrelativistic limit r 2 Foldyfl « r Eh , in QCD both constituent masses and hadronic radii 
are determined by Aqcd so they are expected to be of comparable magnitude. 



III. CONCLUSIONS AND FINAL REMARKS 



The principal conclusion of this paper is that, within the context of the valence quark 
model, the apparent contribution 3/i n /2m^ of the Foldy term to the charge radius of the 
neutron is illusory: in the leading approximation it is exactly cancelled by a "nonintuitive" 
contribution to the radius r\ n of the Dirac form factor FJ 1 . It is therefore totally appropriate 
to compare the measured r 2 En and in general G E (Q 2 ) against quark model predictions (see, 
e.g., Fig. 1 of the second of Refs. [[|) for the rest frame internal charge distribution of the 
neutron. 

Before too much is made of this successful prediction of the valence quark model, some 
other very fundamental questions must still be answered. Perhaps the most fundamental 
is the possible effect of nonvalence components in the neutron wavefunction. After all, the 
classic explanation [pl| for r 2 En is that the neutron has a pn~ component in its wavefunction 



|17| (for a discussion in the more modern context of heavy baryon chiral perturbation theory, 
see Ref. [|18j). Since both hyperfine interactions and qq pairs are 1/N C effects, I know of no 
simple argument for why one should dominate. 

Fortunately, there is both theoretical and experimental progress in resolving this old 
question. Recent theoretical work on "unquenching the quark model" |T9[ indicates that 



there are strong cancellations between the hadronic components of the qq sea which tend to 
make it transparent to photons. These studies provide a natural way of understanding the 



successes of the valence quark model even though the qq sea is very strong, and in particular 
suggest that the precision of the OZI rule is the result of both a factor of 1/N C and strong 
cancellations within this 1/A^-suppressed meson cloud. New data on the contributions 
of ss pairs to the charge and magnetization distribution of the nucleons is also beginning to 
constrain the importance of such effects and, by broken SU{3), their uu and dd counterparts, 
and future experiments will either see ss effects or very tightly limit them (at the level of 
contributions of a few percent to r\ and /zjv). The resolution of the old question of the origin 
of \i n and r 2 En is thus within sight. 

REFERENCES 



[1] R.G. Sachs, Phys. Rev. 126, 2256 (1962). 

[2] R. Carlitz, S.D. Ellis, and R. Savit, Phys. Lett. 64B, 85 (1976); Nathan Isgur, Acta Phys. Pol. 
B8, 1081 (1977). 

[3] Nathan Isgur, Gabriel Karl, and Roman Koniuk, Phys. Rev. Lett. 41, 1269 (1978); 45, 1738(E) 
(1980); Nathan Isgur, Gabriel Karl, and D.W.L. Sprung, Phys. Rev. D23, 163 (1981). 

[4] For some variations on the theme of Refs. see e.g., A. Le Yaouanc, L. Oliver, O. Pene, and 
J.C. Raynal, Phys. Rev. Dll, 680 (1975); D12, 2137 (1975); D13, 1519(E) (1976); D15, 844 
(1976); D18, 1733 (1978); A. Buchmann, E. Hernandez, and K. Yazaki, Phys. Lett. B269, 
35 (1991); Nucl. Phys. A569, 661 (1994). For relativistic quark model variations see, e.g., 
H.J. Weber, Phys. Rev. D49, 3160 (1994); Z. Dziembowski, Phys. Rev. D37, 778 (1988); W. 
Konen and H.J. Weber, Phys. Rev. D41, 2201 (1990); F. Cardarelli, E. Pace, G. Salme, and 
S. Simula, Phys. Lett. B357, 267 (1995); and for the results of a relativized quark model, see 
M. Warns, H. Schroder, W. Pfeil, and H. Rollnik, Z. Phys. C45, 627 (1990). 

[5] See, for example, G. Quemener, C. F. Perdrisat, M. K. Jones, V. Punjabi, and the Hall A 



10 



Collaboration at Jefferson Lab, to be published in the proceedings of the International Nuclear 
Physics Conference, Paris (1998); and V. Punjabi et al, BAPS 43, 1592 (1998). Preliminary 
results from this experiment indicate that G P E (Q 2 ) / G P M (Q 2 ) drops appreciably below unity in 
the GeV 2 region. 

[6] For a recent review see F. Klein, "Recent Experiments on the Electric Form Factor of the 
Neutron", in Proceedings of the 14*^ International Conference on Particles and Nuclei, ed. C. 
Carlson and J. Domingo (World Scientific, Singapore, 1997), p. 121. 

[7] S. Kopecki et al, Phys. Rev. Lett. 74, 2427 (1995). 

[8] R.G. Arnold, C. Carlson, and F. Gross, Phys. Rev. C23, 363 (1981). 

[9] Wigner rotation is the name given to the feature that when a free Dirac particle of mass mu 
with transverse momentum kx is boosted by velocity v along z, its spin flips with an amplitude 
proportional to vkT/mn', see H.J. Melosh, Phys. Rev. D9, 1095 (1974) for the introduction of 
such rotations to the quark model. 

[10] L. Foldy, Phys. Rev. 83, 688 (1951);87, 688 (1952);87, 693 (1952); Rev. Mod. Phys. 30, 471 
(1958). For an example of the impact of Foldy 's paper, see, e.g., the discussion on page 439 of 
the textbook by S. Gasiorowicz, "Elementary Particle Physics" (John Wiley, New York, 1966); 
for a more recent example of a discussion of the Foldy term, see F. Klein in Ref. raj. 

[11] See, e.g., James D. Bjorken and Sidney D. Drell, "Relativistic Quantum Mechanics" (McGraw 
Hill, New York, 1964), p. 52, from which 5r 2 = 3/4m 2 follows; for a glimpse of the early 
literature, see T.D. Newton and E.P. Wigner, Rev. Mod. Phys. 21, 400 (1949). 

[12] See, e.g., F. Gross in "Modern Topics in Electron Scattering" edited by B. Frois and I. Sick 
(World Scientific, New Jersey, 1991), p. 219 for the appearance of r 2 N zwitter in the context of 
relativistic nuclear physics. 

[13] See, e.g., N. Isgur and C. Hayne, Phys. Rev. D25, 1944 (1982) for the appearance of r 2 zwitter 
in the context of the relativized quark model. 



11 



[14] See, e.g., J.D. Bjorken in Results and Perspectives in Particle Physics, Proceedings of the 
A th Rencontre de Physique de la Valle d'Aoste, La Thuile, Italy, 1990, edited by M. Greco, 
(Editions Frontieres, Gif-sur-Yvette, France, 1990) and N. Isgur and M.B. Wise, Phys. Rev. 
D43, 819 (1991) for the appearance of t'q zwitter in the context of heavy quark symmetry, 
where it corresponds to the fixed minimum slope of \ for the Isgur- Wise function. 

[15] For the proton J2i e i r i = 3/4mg, which is much larger than the Foldy term and, by contributing 
~ 0.27 fm 2 to r\ p , brings the constituent quark model prediction for this quantity into better 



agreement with experiment [ 13 1 



[16] In its early manifestations, the pion cloud model was posited to give rise to the nucleon form 
factors by pion vertex corrections analogous to the QED diagrams responsible for g— 2. See, e.g., 
S.D. Drell and F. Zachariasen, "Electromagnetic Structure of Nucleons" (Oxford University 
Press, New York, 1961); R. Hofstadter, "Nucleon and Nuclear Structure" (W.A. Benjamin, 
Inc., 1963); L. Hand, D.G. Miller, and R. Wilson, Rev. Mod. Phys. 35, 335 (1963); S.D. Drell, 
in Int. School of Physics "Enrico Fermi", Course XXVI, Varenna, 1962 (Academic Press, New 
York, 1964). As explained in this work, the strong irir interactions in the t-channel call for 
dispersion-theoretic approach. See G.F. Chew, R. Karplus, S. Gasiorowicz, and F. Zachariasen, 
Phys. Rev. 110, 265 (1958), and P. Federbush, M.L. Goldberger, and S.B. Treiman, Phys. Rev. 
112, 642 (1958). These papers in turn led to the development of the vector meson dominance 
model. For a sense of the degree of dominance of the meson cloud picture of the nucleon form 
factors before the discovery of the quarks, see the figure on page 451 of the undergraduate 
textbook "Mechanics", Berkeley Physics Course - Volume 1, by Charles Kittel, Walter D. 
Knight, and Marvin A. Ruderman (McGraw Hill, New York, 1965): the point-like proton 
core and pion cloud of both the proton and neutron are evident and stated to be based on 
high-energy electron scattering data. 

[17] I have examined the Foldy term in a model of a nonrelativistic PD L = 1 bound state, where 
P is a pseudoscalar particle. In such a system (which represents the most naive version of a 



12 



pir~ model for the neutron), p sp in-flip makes a contribution which exactly cancels r 2 D zwitter l 
Thus in the pseudo-charge-conjugation symmetry limit mp = mo = m, r 2 EhpD = and, as in 
the quark model, one has r\n PD = ~ r Foidy h PD so ^ na ^ an y observed r 2 En may be interpreted as 
arising from an intrinsic internal charge distribution (as naively pictured in the pir~ model) . 

[18] See, e.g., V. Bernard, H.W. Fearing, T.R. Hemmert, and U.-G. Meissner, Nucl. Phys. A635, 
121 (1998) and references therein. 

[19] P. Geiger and N. Isgur, Phys. Rev. D 41, 1595 (1990); P. Geiger and N. Isgur, Phys. Rev. D 
44, 799 (1991); Phys. Rev. Lett. 67, 1066 (1991); Phys. Rev. D 47, 5050 (1993); P. Geiger, 
ibid. 49, 6003 (1993); P. Geiger and N. Isgur, Phys. Rev. D55, 299 (1997). 

[20] "Measurement of the Neutral Weak Form Factors of the Proton" , The HAPPEX Collaboration, 
K.A. Aniol et al., submitted to Phys. Rev. Lett. 



13 



